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' In this paper we study the production of energetic di-leptons. We calculate the rate for 2 — > 2 

' processes. The log term is obtained analytically and the constant term is calculated numerically. 

, When the photon mass is of the order of the thermal quark mass, the result is insensitive to the 

^ ■ photon mass and the soft logarithmic divergence is regulated by the thermal quark mass, exactly as 

Q ' in the case of real photons. We also consider the production of thermal Drell-Yan dileptons (thermal 

quark and antiquark pairs produced by virtual photons) and calculate the rate systematically in the 
context of the hard thermal loop effective theory. We obtain analytic and numerical results. We 
' compare our results with those of previous calculations. 
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I. INTRODUCTION 



It is believed that experiments at RHIC are producing quark-gluon plasma. The properties of the plasma will be 
further studied by the heavy ion program at the LHC. Thermalization and rescattering tend to erase information 
about the state of the plasma at early times, and therefore it is of interest to study high energy photons which are 
produced early in the collision and escape without further interaction (for a review, see jl|). 

The main disadvantage of using photons as a signal is that there are several sources of photons, and it is necessary to 
separate the signal from the background. Photons that are produced in the early stages of the collision are called direct 
photons. There are two types of direct photons: prompt photons, which are produced in the initial collisions of the 
partons which make up the heavy ions, and thermal photons, which are produced in the hot quark or hadronic matter 
' _ formed during the collision. Prompt photons have a power damped spectrum and dominate over thermal photons at 
I ■ large energies. Thermal photons produce an exponentially damped spectrum and dominate prompt photons at lower 
— . ' energies. In addition to direct photons, there is a large background contribution from photons produced through the 
. radiative decay of hadrons. These photons do not provide direct information on the early stages of the quark-gluon 
' plasma. 

• ^ , Of these three kinds of photons (prompt /thermal/decay), we are primarily interested in thermal photons, since they 
' are the ones that provide direct information on the early stages of the plasma. In order to isolate the contribution 
H I of thermal photons, one should look at not too large transverse momentum, typically a few GeV's. To overcome the 
- - < background problem in this range, it is of interest to study the spectrum of small mass di-leptons (virtual photons) 
which have the same production mechanisms as real photons. In fact, the PHENIX experiment at RHIC has been 
able to extend the accessible qt range of the photon spectrum down to 1 GeV using this approach Q . Throughout 
this paper we use go for the photon energy, Q := {Q^ = Qa ~ 1^) for the invariant mass, and T for temperature. 

We focus on the case qo ^ T ^ Q. 

In general, the production rate for real or virtual photons is obtained from the imaginary part of the contracted 
retarded photon polarization tensor (see Eqns. (HI) and ([2])). We work in the framework of the "Hard Thermal Loop" 
(HTL) resummed perturbation theory Q where the quarks and gluons acquire an effective mass as well as a dispersive 
component. 

For real photons, there are two physically distinct contributions to the polarization tensor at leading order. 

Real photons are produced through Compton scattering and 2 — > 2 annihilation processes. Collectively, these will 
be called 2^2 processes. They are obtained from the central, or real, cuts of the two loop diagrams, where the gluon 
is time-like and on shell. The 2-loop contributions to the polarization tensor are shown in Fig. [T] The amplitudes 
that are obtained from the real cuts with time-like gluons are shown in Fig. [51 



FIG. 1: 2-loop contributions to the polarization tensor; solid lines are quarks, dashed lines are photons and dotted lines are 
gluons. 



FIG. 2: Amplitudes that correspond to 2 2 processes. Crossed diagrams are not shown. 

In addition, there is a leading order contribution to the rate from bremsstrahlung and off-shell annihilation processes. 
These processes contribute at leading order because of a strong coUinear enhancement. They arc obtained from the 
real cuts of the two loop diagrams in Fig. [TJ where the gluon is taken to be HTL resummed and space-like. The 
scattering amplitudes are shown in Fig. [31 




FIG. 3: Amplitudes for bremsstrahlung and off-shell annihilation. 

For the case of space-like gluons, higher order diagrams involve the same enhancement and therefore also contribute 
at leading order. This is known as the LPM effect. The diagrams are resummed using a set of integral equations (^j- 
Real cuts with space-like gluons produce a set of amplitudes that corresponding to multiple rescatterings. For off-shell 
annihilation processes, some of these amplitudes are shown in Fig. 31 The rate obtained from Eqn. ^ involves terms 
that come from squaring these amplitudes, and also terms that correspond to interference effects. 
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FIG. 4: Some multiple rescattering diagrams. 

Now we discuss the case of virtual photons. In this case a new scale Q enters the calculation. The complete 
calculation has already been done, but due to several misunderstandings, the pieces have not been correctly assembled. 

We start with the observation that in HTL the quark acquires an effective mass ttt-oo , and the imaginary part of the 
1-loop diagram has a threshold at > 4to^ for a thermal Drell-Yan processes of the form qq~^ ^, which corresponds 
to a thermal quark anti-quark pair produced from a virtual photon. The self energy diagram and the amplitude that 
is extracted from the imaginary part are shown in Fig. [5l [l5|. 




FIG. 5: The lowest order thermal Drell-Yan process. 
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There will also be contributions from the side cuts of the 2-loop diagrams shown in Fig. [1] and from the side cuts of 
all contributions to the polarization tensor that produce the amplitudes shown in Fig. [D When the gluons are taken 
to be space-like, these side cuts give contributions that correspond to interference between the tree level amplitude 
and virtual corrections to it. Some of the amplitudes are shown in Fig. [S] (again, we show only annihilation processes). 
The dilepton rate obtained from Eqn. ([T]) involves terms that come from squaring these amplitudes, and also terms 
that correspond to interference effects. 
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FIG. 6: The amplitudes involved in the LPM resummation for virtual photons 



In principle, there is also a contribution at leading order from the side cut of the 2-loop diagrams where the gluon 
is time-like. The contribution from soft time-like gluons is subleading. The contribution from hard time-like gluons 
is included by using asymptotic propagators and vertices in the self energy diagrams that produce the amplitudes 
shown in Fig. [6l This point will be discussed in detail in section Hvl 



The results for these contributions that have been obtained by previous calculations are summarized below. 
For real photons, the rate for 2 2 processes has been calculated in Refs. [1,0]. 

For real photons the rate for bremsstrahlung and off-shell annihilation processes has been calculated in Ref. [1]. 

For virtual photons the rate for bremsstrahlung and off-shell annihilation processes has been calculated in Ref. 0] ■ 
The leadingterm in the LPM resummation is the rate for the thermal Drell-Yan ^ 7 process. 

In Ref. [8[ the thermal Drell-Yan contribution was calculated to logarithmic accuracy assuming go ^ Q S> to/ 
(where m/ is the quark thermal mass). A comparison with our calculation is given in section IVAI 

The thermal Drell-Yan contribution was also calculated in Ref. [9|. These authors attempt to produce a result 
that is valid at next-to-leading order for qq ^ T ^ {m/,Q}. There are several problems with this calculation. A 
detailed discussion is given in section IVBl 

This paper is organized as follows: In section we define our notation. 

In section Hill we calculate the rate for 2^2 processes for virtual photons. We include the new scale Q, and an 
asymptotic mass for the hard fermion. We show that the asymptotic mass does not affect the result. The log term is 
obtained analytically and the constant term is calculated numerically. Our result correctly reduces to the result for 
real photons, in the limit that the photon mass goes to zero. 

In section IIVI we calculate the contribution to the rate from thermal Drell-Yan processes both analytically and 
numerically. The numerical result is valid for qq ^ T ^ {'^/i Q}j without any restrictions on the relative size of Q 
and ruf. We produce two analytic expressions: one is valid for small mf and one is valid close to the threshold. 

In section IVl we compare our results with those of previous calculations. 

In section IVll we present our conclusions. 



II. NOTATION 



In this section we define our notation. We write Imll := Im(nrct)^^- The rate for virtual photons (di-leptons) is: 



rate :— — — — = 7— rlmll. (1) 

d^xd^q 127r4Q2 ^ 



The rate for real photons is given by a similar expression: 



d^T e~'^'^'> 

■= ■ d^q = (2^^"^n- ^2) 



We use conventional notation for the thermal distribution functions: 
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We write 4- vectors as capital letters, for example P = {po,p). In addition, we give specific names to particular 
combinations of momenta (see Fig. [1]): 



S = L + P + Q, T = L + P, R = P + Q, S' = -S. 



We define angles as follows: 



q= (0,0, g), 

p—{p sin /3, 0, p cos /3) , 

s' = (s' cos (/()sin7, s' sin (/)sin7, s' COS7). 

We write the cosines of these angles as: 

COS P = Xjs , cos J = x^. 
For bare on shell fermions and gauge bosons (in the Feynman gauge) we write: 

^rct(P)-5adv(P) =f d(P), 

DZ{L)-D^^^^iL) = -g^-^diL), 
d{X) = -27riSign(2;o) S{X^). 

For a HTL fermion: 

SrotiP) - 1 ((70 -pV)A+t(^) + (70 +pV)Arct(^)) 

The spectral density is obtained from: 

P±{P) = *rf±(P) = *(A+t(P) - A+ ,(P)), 



(4) 



(5) 
(6) 



(7) 



(8) 



d{P) 
Z± 

f3+{po,p) 



-2m ^ Z±S{pa - juj±) - 2TTif3{po,p), 



(jjj- — p 



2p^m) {p ~ Po) Q {p^ - pI) 



{p - Pof mj + (2{p- po) p2 + mj (^p 



In 



P+Po 
po-p 



2 - In 



P+Po 
po-p 



Po 



III. COMPTON AND ANNIHILATION PROCESSES 



We expect that leading order contributions to Compton and annihilation processes come from the real, or central, 
cuts of the diagrams in Fig. [U where the gluons are time-like. The first diagram in this figure is referred to as the 
self-energy diagram and the second is called the vertex diagram. We note that there are two self-energy diagrams, 
which correspond to correcting either the top or bottom fermion line. Both diagrams give the same contribution and 
thus we could use either one of them with a factor of two but, in order to obtain a more symmetric expression, we 
explicitly include both diagrams. 

Only the hard part of the gluon momentum contributes, and thus it is safe to use a bare on-shell gluon propagator. 
On the other hand, both hard and soft quark momenta contribute. In order to calculate the integral, we separate the 
phase space into two regions by introducing a scale Pc ~ ^/gT. We calculate separately the contributions involving 
exchanged quarks with momenta greater than and less than pc- Different approximations are used in each region. 
The separation scale Pc cancels when the two pieces are combined. This technique was introduced in [lo| . 

At zero temperature, the integral corresponding to the diagrams in Fig. [T]can be written: 



Imn = 



(27r)8 



dP / dS 



(9) 



(2Tr(5(P)7p5(i?)7.5(5)7A5(i?)7'0+Tr(5(P)7.5(r)7^5(5)7A5(i?)7^))i^^'(i), 



5 



where = 5/9, Cp = 4/3, Nc = 3, J dP := J dpo f d^p, and iS and —iD^i, correspond to the quark and gluon 
hnes respectively. We need to obtain the corresponding integral at finite temperature. We work in the Keldysh 
representation of the real time formalism. The method we use to sum over Keldysh indices is described in [llj and 
our technique for taking the imaginary part is found in . We calculate the contribution for hard exchanged quarks 
by explicitly imposing a lower limit cutofl^ on the quark momentum. The calculation is completely insensitive to soft 
momentum scales, which means that we can use bare propagators (Eqn. ([6])). In addition, we drop terms in the 
numerator that are proportional to , since we have assumed go 3> {Q, rnj} throughout. The result for the real cuts 



is: 



Imn = " ^ ' ^''^ J dpp^ J dpo J ds' {s'f j ds'o J dxp j dx^ J d0 [SE + VER] (10) 
• dLds'dp ((1 + Ub {lo)) (1 - n-/(so)) (1 ~ (Po)) - nt {lo) nf{sQ)nf (po)) , 



2V 



SE = ^ -r—, '- ■ VER = 



i?2r2 ' i?2y2 

The term SE corresponds to the first diagram in Fig. [1] combined with the diagram with the bottom fermion 
corrected. The term VER corresponds to the second diagram in Fig. [T] The result is exactly the same for real and 
virtual photons. Following Q we obtain: 



T^rrhard o AT 2^2 /i / \ , ^hard\ /-fhard C'(2) ln(2) 3 

Imn2^2 = 2aiVce,my lln I — 1 +0*2^2 I ; '^2^2=-^^ + —^ ^+2' 

For soft exchanged quarks a slightly different technique is necessary since the integral diverges when the virtuality 
of the exchanged quark [R^ or T"^) goes to zero. This is just the expected result that the cross section involving 
the exchange of a massless particle is infinite. The well known solution to this problem is to replace the soft quark 
propagator with the corresponding HTL propagator. Effectively, when the exchanged quark is soft, the diagrams in 
Fig. [T] are replaced with the diagrams shown in Fig. [71 where the propagator with the solid dot is the HTL quark 
propagator and the bare line is a hard quark propagator. Using symmetry, we calculate the first diagram and multiply 
by a factor of two. 



FIG. 7: 1-loop contributions to Compton and annihilation processes. The quark lines carry momenta P and R = P + Q. 

The integral corresponding to two times the first diagram in Fig. [7|is given in Eqn. (|12p . The dressed line carries 
momentum P and the hard line has momentum R = P + Q. 



Imn = 



4z7raiVf,e2 



dpo / dpp'^ / dxfj {nf (po) ~ Hf {po + qo)) 5 (^{po + qof ~ i^irf'j (12) 



id+ {po,p) {-p + Po + qo ~ qxfj) +d_ (po,p) {p + Po + qa +qxp)) 

The real cut is obtained from the cut part of the HTL dispersion relation d±(jJo,p) — 27ri/3±(po,p). Since there is 
an explicit factor rrij ~ ctsT^ in the numerators of the factors P±(j>o,p) (see Eqn. ([8])), we will drop all additional 
factors of order {m/, Q} in the numerator of (fT2|) . 

We use the delta function to do the angular integral. Setting the argument of the delta function to zero and solving 
for Xi3 gives the solution Xfj. Writing pq = zp we define g{z,p) := 1 — x'^ and obtain the phase space constraint 
g{z,p) > 0, which we can write as a theta function. The limits on the p-integral are determined from this constraint. 
We use u}{r) = ^/r"^ + m'^ = •y/p2 + 2qpx/3 + q'^ + rri^ where we define the asymptotic quark mass by = 2to2^ 
Doing the angular integral we obtain: 

Imff2°l*2 =8aN,elm} J^' J ^dzdp {nf{pz)-nf{pz + q„)) {f{z,p) + f{-z,p))e{g{p,z)), (13) 

p^{z + lf 



f{z,p) 



7r2(z + l)2m4. + (^{z + l)ln (^f±i^ m2 - 2 (^{z + l)p^ + mj^y 
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where we have defined: 



Since T ^ pc > p we can use 



g{z,p) = - - 1) p2 + 2zqop + - 



nf{pz) - nf{pz + qo) 



nfiqo)) 



(14) 



(15) 



to factor the thermal functions out of the integrand and obtain: 



dzdp {f{z,p) + f{-z,p)) e{g{p,z)) 



(16) 



Corrections to this result are or order Pc/T and are thus of the same order as terms that have been dropped in the 
calculation of the contribution from hard exchanged quarks. Recall that we are working in the limit go ^ S> Pc 3> 
{uifjQ}^ and it is only in this limit that the dependence on the cutoff cancels between the contributions from hard 
and soft exchanged quarks. 

Using go ^ {"^/i Q} we can rewrite the constraint given by the theta function in terms of integration limits: 



490 



<P<Pc 



and 



-1 < z < 1 - 



2qop 



for g2 - > 



(17) 



2<?oP 



< z <1 for TO^ - > 



Using these results, for fUf/qo <^ 1, the integral is completely parametrized by the dimcnsionless variable ^ : = 
gV(gom/). 

The logarithmic term is extracted by taking the static limit z ^ 0. The plus and minus modes give the same 
contribution to the log term. The result is: 



Imn^'^*2 = 2aNcel ml 



In 



^0 J 



(18) 



where Cl%{Cj is a constant that must be determined numerically. Combining (fTT|) and we see that the cutoff 
Pc cancels. We obtain: 



ImUs, 



2 a TVc mj 



■In , 

2 



Chard i ^^soft 
9 .0 I O9 ,9 



(19) 



We consider the denominator of the log in Eqn. (I19|) . For gn ^ Q ~ ti/ (^ ^ 1) we can drop the second term in the 
square brackets and the expression reduces to the result of 0,0] for real photons. However, we must be careful when 
^ is large. It is easy to see that the terms that are proportional to the asymptotic mass are always negligable (by 
powers of Q/qo) compared to terms proportional to the thermal mass. This result indicates that it is not necessary 
to include the asymptotic mass on the hard line in Fig. [71 or equivalently, that the coUinear part of the integral is 
not dominant. Dropping the asymptotic mass we obtain: 



Inin2^2 = 2aNcelmj 



-In. 

2 \^mj 



Note that for ^ 1 all three terms in the denominator of the log are of the same order. 



(20) 



For ^ = the numerical result for the constant C|2l*2(C) agrees with the result for real photons @, [1]. 

For ^ ^ 1 we can obtain an analytic expression (up to corrections of order Q^/goPc) for the integral in Eqn. (fTC)) : 



ImU? 



= 2 a TVc mj 



dp 

QV(4go) P 
In 



1-QV(2P90) 



dz 



(21) 



Aplql 



2 In 2 - 2 
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Comparing Eqns. ^ and ^ we obtain C^'^^{^ ^ oo) = 21n2 - 2. 
The full numerical result for the constant C|"£*2(0 shown in Fig. [H 




FIG. 8: The term C|!il2(0 as defined in Eqn. ([18}. For ^ ^ oo the result approaches 21n2 - 2 ~ -0.61. For ^ = the result is 
approximately -0.69. 



IV. THERMAL DRELL-YAN PROCESSES 



A. Integral Expressions 

In this section we discuss the contribution from thermal Drell-Yan processes of the form qq 7. For these 
processes we have a hard, low virtuality photon {{qo, q} hard; Q soft) decaying into two quarks with momenta P and 
R ~ P -\- Q. By momentum conservation, at least one of these quarks must be hard. In fact, because the virtualities 
of both quarks are required to be soft by the kinematics, the integral is phase space suppressed unless both quarks 
are hard {{pQ,p,rQ,r} hard; {P^,i?^} soft). The behaviour of hard, low virtuality quarks was first studied in p^ . 
Since we are effectively extracting a next-to-leading-order contribution from a one-loop diagram, we need to keep all 
next-to-leading-order terms in the asymptotic propagator. Following Ref . , we write: 

5'htl(^) -> 5'asy(-R) = p2 ' ^^^^ 



It is normally sufficient to use bare vertices when all legs are hard. However, when calculating a next-to-leading- 
order contribution, one must also use corrected vertices. One way to understand this point is to notice that the 
asymptotic propagator and the bare vertex do not satisfy the Ward identity. 

Combining this information, we consider the 1-loop diagram in Fig. [5) 




FIG. 9: The diagram from which thermal Drell-Yan processes are obtained. The quarks carry momenta P and R = P + Q. 
We use the asymptotic propagators given in Eqn. (|22|l . and the vertices given in Eqn. (|23|l . 



8 



We use corrected vertices of the form [T3 |: 

= {g^.u + ST'f^ir (23) 

The exphcit form of the vertex correction is not necessary. We only need to require that it satisfies the Ward identity: 

ST'^l ^P,-P,- 5r;;t' ^Ru-Ru- (24) 



The integral corresponding to the diagram in Fig. [9] is constructed from Eqns. ([22|) and ((24|) . 

We remark that the rate for 2 ^ 2 processes is insensitive to asymptotic masses and vertex corrections because the 
kinematics does not require both of the quarks with momenta R and P to be on-shell. We have verified this by direct 
calculation. 

The imaginary part of 11 is obtained by cutting both lines. We do the angular integral using the delta function 
5{{P + QY' — rn^) as in section ITTll After using the second delta function to do the po integral, the argument of the 
theta function in Eqn. (|14p becomes: 



g{-tOoo{p)/p,p) - 4p2g2 _ (q2 _ 2go^oo(p))' (25) 
and determines the limits Pmax and Pmin on the p-integral: 

Pmax = ^ (? + goVl - 4to^/Q2^ , Pmin ^ 9 - ^0 V 1 - 4m^/Q2 . (26) 

The condition Pmax — Pmin immediately shows that the threshold for the production of a thermal Drell-Yan pair is 

The integral corresponding to the diagram in Fig. [5] is given by: 

2aelNc f 

ImHoY = 1— dP {nf{po)~nf{po + qo)) (27) 



• d{P^ - m^) d((P + Qf - m^) [P-{P-P)+R-{R-R)-P-R] 

where we have dropped terms of order and higher in the numerator, since they do not contribute to the rate at 
next-to-leading-order. Doing the po and angular integrals and writing the result in terms of the variable uj := lUoo(j>) 
we obtain 

ae^,Nc / o mj {Q^ + 2uj {uj - qo)) (u;^ + [qo - ujf) _ 
IrnHDY =^^/ div {nf{~Lj)-nf{qo^u;))\Q^ + — \^ ^ (28) 



^2 (go _ u;) 



where 



i (go + 9v/l-4m^/Q') ; ^min = ^ (90 - - 4rr74/Q') • (29) 



B. Analytic Expressions 

In this section we derive some analytic approximations to Eqn. (j28p . We are interested in two limiting cases: 
qo ^ Q ^ ruf and qo^ Q ^ nif. The first limit will give the behaviour of the thermal Drell-Yan contribution close 
to the threshold. The second limit allows us to compare with previous calculations [1, Q. 



1. Threshold Expansion 

For qf) ^ Q ^ vrif we have tj^max ~ ^min ~ 'Zo/2 and thus cu is hard throughout the full range of the integration, 
which is consistent with our use of asymptotic propagators in Eqn. (|27p . We can expand ((28)) in {qo,^} ^ {f^fjQ} 
and approximate 

Hf {-Ll>) - Uf (qo - uj) sa 1. (30) 
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We obtain: 



duj g2 - 



uj {qo - uj) 



(31) 



This result is the same as Eqn. (20) of Ref. Q This expHcitly demonstrates that the contribution calculated 

in this section (and shown in Fig. [9|) corresponds to the first term in the LPM resummation shown schematically in 
Fig. [51 where the vertices and quark propagators include asymptotic corrections. 
The integral in Eqn. (|3ip can be done analytically. Using Eqn. ([29]) we obtain: 



ImHoY = ae^N, ^(1 " 4m^/Q') (2m^ + Q^) + 2m^ In 



1 + y/l - 4mg,/Q^ 
1-^1- Aml,/Q^ 



(32) 



In order to discuss this result, and to compare with the result of other calculations, we separate the Born term from 
the complete thermal Drell-Yan contribution. We define: 



ImnBorn « aNcel^yq^ {Q^ ^ 4ml,). 



and write the correction to the Born term as: 



ImllpY := ImllDY - ImllBorn- 



We obtain: 



Imnjj^l = iaN^elmj ( y/l - Am^/Q^ - In 



1 + y/l - 4mg,/Q^ 
1-^1- 4ml/Q^ 



(33) 
(34) 

(35) 



2. Small nif Expansion 



In the limit qo ^ Q ^ mf we have Wmax ~ 9o ! Wmin ~ Q^/(4(7o)i which corresponds to Wmax hard and Wmin 
soft. The integral expression corresponding to the thermal Drell-Yan contribution (Eqn. psp ) was derived under the 
assumption that the leading order contribution comes from w-hard. As a consequence, we cannot get a result beyond 
logarithmic accuracy in the small mf limit. Since the entire thermal Drell-Yan contribution goes to zero as my — > 0, 
logarithmic accuracy is sufRcient. We study this limit mainly to compare with the results of other authors. 

In order to extract the log term in the small mf limit we can expand Eqn. (j28p in qo 3> {i^, mj, Q}. Dropping the 
Born term we obtain: 



imn 



(1) 

DY 



mf — >0 



= -2aeiNc 

q 1- 



diO — 



q^ 



(36) 



C. Numerical Results for the Thermal Drell-Yan Contribution 



In Fig. [TOlwe plot the results for the thermal Drell-Yan contribution to Imll. We use 

qo = 5T, T = l GeV, Q = 0.25 GeV. (37) 

In order to clarify the difference between the different calculations, we plot the correction to the Born rate only. 

The solid black line is the exact numerical result for the thermal Drell-Yan contribution Il^y obtained from Eqn. 
([28]) by subtracting the term that gives the Born contribution. The red dashed line is the result of patching together 
the analytic expressions in Eqns. ([35]) and ([36]) . The blue dot-dashed line is the result of Ref. Q. The black dotted 
line (the lowest line) is the result of Ref. (8|. In sections IV Al and I V Bl we discuss in detail the difference between these 
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results and ours. 



Full numeric 
Analytic approximation 
Thoma & Traxler 
Altherr & Ruuskanen 



-0.0015 L 



0.5 
2m„ 



FIG. 10: 



(1) 



0, where II^y goes to zero. Our result 



has the correct threshold at Q = {2m, 



V. COMPARISONS WITH OTHER CALCULATIONS 

In this section we compare our results with the results from some other papers. This task is made more difficult by 
the fact that different notation is used everywhere. A comparison of the definitions of 11 and the rate in the relevant 
papers is given in Eqn. (1381) . The subscript us refers to the definitions used in this paper, refers to the work of 
Aurenche at al (Refs. (tI. I14|)- ^i? refers to the work of Altherr and Ruuskanen (Ref. [8]); and TT refers to the work 
of Traxler and Thoma (Ref. 



dileptons : raters = rateAGZ = 1/2 • rate ar — 1/2 ■ ratcTT, 
Imll : Imllus = -Imll^cz = Imll^ij = ImllTT- 



(38) 



Throughout the rest of this section, results from other papers are translated into our notation. 

A. Comparison with the work of Altherr et al. 

In section Hill we use an HTL propagator for the soft exchanged quark to calculate the contribution from 2 — > 2 
processes, and in section IIVI we use asymptotic propagators to calculate the contribution from thermal Drell-Yan 



processes. The authors of 8] use expanded versions of these propagators by replacing ^ 
consider only massless quarks. Writing only the log terms, their results are: 



imn 



(1) 

DY 



2 a iVc mj In 



ImUo 



^ + and 



(39) 



where the factor m^cg is an arbitrary, unphysical regulator on the quark virtuality. The cancellation of this regulator 
taken to be is evidence of the KLN theorem. The results in Eqn. ((551) ^^'S equivalent to the mj — > limit of our 
results (see Eqns. (PUI and ([551) 1. if the regulator taken to be rrircg = Q^/qo- 



B. Comparison with the work of Thoma et al. 

The authors of intend to produce a result that is valid for ^ T ^ {mf,Q}. They calculate the thermal 
Drell-Yan contribution from a 1-loop graph, with one hard bare line, and one HTL propagator which is taken to be 
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asymptotically hard. The diagram is multiplied by a factor of two, to account for the fact that either propagator 
could be the hard bare one. There are several problems with this procedure. 

Failure to use the asymptotic propagator on both lines gives the wrong symmetry factor. As a result, the Born rate 
(Eqn. (10) in Ref. [9*]) is too big by a factor of two. In the result for the correction to the Born rate, the factor of 2 
was inadvertently dropped. 

In addition, using one hard bare propagator gives the wrong threshold. The limits on the p-integral are obtained 
using a theta function like the one we have in Eqn. (j25p . Using our notation, the argument of their theta function is: 

- (/ - 2qp - Q2 _ ^^(p)2 + 2qoLO+{p)) (/ + 2qp - - uj+{pf + 2qoLO+{p)) (40) 
Analytic approximations are obtained using the asymptotic dispersion relation uj+ip) = p + 'm'j/p which gives: 

2m? go 

Pmax ~ qo ; Pmin ~ ~~Q2~ ' ^^^^ 

From these expressions, one can see immediately that the threshold obtained from Pmax — Pmin occurs at = 
^rnj = m^. One can also find the threshold numerically using the full HTL dispersion replation in (PD|) . The result 
is virtually unchanged. This threshold corresponds to a decay into one quark with mass rriao and one massless quark, 
or the annihilation of a massive and massless quark pair. The correct result is = (2moo)^, which is what we get 
by including the asymptotic mass on the hard line. 

Finally, by neglecting vertex corrections the authors have missed some leading order contributions. 



C. Discussion of the results of Aurenche et al. 



The full virtual photon spectrum has the following components: 

(1) The 2^2 processes as given in Eqn. ([^D]) . We write this contribution schematically as M2_>2 (see Fig. [2]). 

(2) The thermal Drell-Yan term and the associated off-shell annihilation and bremsstrahlung processes, including 
multiple scattering contributions and interference terms, all of which are contained in the LPM resummation which 
was calculated in We write schematically the matrix elements corresponding to these contributions (see Fig. [5]): 

oo oo 
n=0 n=0 

The first term in the annihilation sum is A/2^'y and corresponds to the contribution from thermal Drell-Yan processes. 

Refs. 0, dealt with bremsstrahlung and off-shell annihilation processes for virtual photons (see Fig. [S]). In 
order to compare the contributions to the overall rate from various terms, an estimate of the 2 —^ 2 processes was 
constructed in the following way. The contribution from 2^2 processes was taken to be the result obtained by 
Altherr and Ruuskanen [8j. This result contains a \n{qQT/Q^) term and thus for small enough the result of \^ 
becomes larger than the rate calculated in [5, 6] for real photons (which contains a logarithmic factor ln{qoT / m'j)) . 
It is clear that this is unphysical: when Q ^ we drop below the threshold for thermal Drell-Yan processes and the 
rate for virtual photons should reduce to the result for real photons. Therefore, for small Q^, the result of [H was 
replaced by the result of d, 0] for real photons. As a result of this patching procedure, the rate as a function of 
contains a slight corner at the value of where the two curves cross (see, e.g.. Figs. 6 and 7 of Ref. 14]). Clearly 
this corner is an unphysical feature. 

We note that the rate calculated by Altherr and Ruuskanen is the small to/ limit of the properly calculated 
thermal Drell-Yan process with resummed propagators and associated vertex corrections (Eqn. (j28p ). Furthermore, 
as discussed above, the dressed thermal Drell-Yan process is the first (order 0) term in the integral equation which 
resums the re-scattering corrections for bremsstrahlung and off-shell annihilation, and is therefore already included 
in the calculation of 0] . The thermal corrections to the Drell-Yan process have therefore been inadvertently included 
twice in ^, ^ij] . 

The correct expression for the rate is the sum of the result for the LPM resummation which was calculated in [3], 
and result for 2 ^ 2 processes as given in Eqn. ([^0)1 . 



VI. CONCLUSIONS 



We have calculated the contribution from 2 — > 2 processes to the production of virtual photons (Eqn. (PO)) ). Our 
result is valid for arbitrary values of the parameters qo 3> {to/, Q}. The log term is obtained analytically and agrees 
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with the result for real photons The constant term is obtained numerically and shown in Fig. [S] In the limit 

Q ^ it agrees with the result for real photons. 

We have also calculated the contribution to the rate from thermal Drell-Yan processes both analytically and 
numerically. The numerical result is valid for qo ^ T {mf,Q}, without any restrictions on the relative size 
on Q and mj. Wc produce two analytic expressions: one is valid for qo ^ Q ^ mf and one is valid close to the 
threshold. 

The thermal Drell-Yan rate has been calculated twice previously (Refs. ^1]), but neither calculation is completely 
correct. The thermal Drell-Yan contribution was also calculated in Ref. 1^, as the leading order term in the LPM 



resummation of bremsstrahlung and off-shell annihilation processes. A slight error occured in the presentation of 
these results, due to a misinterpretation of the results of [sl,!^]. 
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The threshold for processes of the form q q + ^ is < (mi —mT)''. As is discussed in section II VI both quarks carry HTL 
propagators where the momentum is taken to be much larger than the thermal quark mass. In this limit, iu+ « \/jP'^\-rn^ 
and oj- ~ p and therefore the process q ^ q + "/ must involve one plus mode and one minus mode. Since the residue of 
the minus mode is exponentially suppressed at large momentum, there is no phase space for this process. 
[16] The authors of [Tj] define the polarization tensor with the opposite sign. Also, there is a missing factor l/(2g) in their Eqn. 
(20). 



